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Abstract 

We strengthen and generalise a result of Kirsch and Simon on the 
behaviour of the function Nl(E), the number of bound states of the 
operator L = A + V in K d below —E. Here V is a bounded potential 
behaving asymptotically like P(u>)r~ 2 where P is a function on the 
sphere. It is well known that the eigenvalues of such an operator are 
all nonpositive, and accumulate only at 0. If the operator A s< j-i + P 
on the sphere has negative eigenvalues —fii, ■ ■ ■ , less than — {d — 
2) 2 /4, we prove that Nl(E) may be estimated as 



N L (E)) = £ " (d ~ 2) 2 /4 + O(l); 



thus, in particular, if there are no such negative eigenvalues then L has 
a finite discrete spectrum. Moreover, under some additional assump- 
tions including that d = 3 and that there is exactly one eigenvalue 
—fj,\ less than —1/4, with all others > —1/4, we show that the neg- 
ative spectrum is asymptotic to a geometric progression with ratio 

exp(-2vr/ J m - \). 
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1 Introduction 



Consider a Schrodinger operator of the form H = A + V acting on L 2 (M. d ) for 
d > 3, where A = — ^ df is the positive Laplacian and V is a multiplication 
by a potential function V(x) which tends to zero at infinity. It is a basic 
problem in quantum mechanics to determine the location and the nature 
of the spectrum of H. It is well known, and not hard to show, that on 
(—00, 0) the spectrum is discrete. It is known that under rather general 
conditions (decay of V at some prescribed rate at infinity, (^(r" 1 ) being 
sufficient), there are no positive eigenvalues of H jH], [0]. Finally, under 
fairly general conditions, e.g. V = 0(r _1_e ) for some e > 0, the restriction 
of H to its continuous spectral subspace is unitarily equivalent to A (see e.g. 
Kato |H] or Yafaev 15J, Theorem 2.4). All of these results apply under the 
assumptions on V made below. The main question, related to the description 
of the spectrum at this coarse level, left unanswered by the considerations 
above is whether the point spectrum accumulates at 0, i.e. whether the point 
spectrum is finite or infinite. Finding bounds or asymptotics on the negative 
spectrum has been a major topic of research in Schrodinger operators for 
many decades, and it remains an active field (see e.g. [Hj, [Z]). 

Let us denote by Nh(E) the number of eigenvalues of H below —E < 
(assuming inf a ess (H) =0). It is well known that the behaviour of V near 
infinity determines whether Nh{E) is bounded as E — > 0. There is an 
attractive discussion of this in jT3|, section XIII. 3. For example, if V > 
_ cr -2-e f Qr some positive c then Nh(E) is bounded, and it is unbounded 
if V < —cr~ 2+e (here r = \x\ is distance to the origin). In the borderline 
case where V behaves asymptotically as — cr~ 2 , then Nh{E) is bounded if 
c < l/4(d - 2) 2 (in fact, the operator A — (d — 2) 2 r" 2 /4 on L 2 (R d ) is a 
positive operator) and unbounded if c > (d — 2) 2 /4, where d > 2 is the 
spatial dimension [TUj . 

In this paper we shall analyze the more general case where V ~ cr~ 2 along 
any ray from the origin, but the constant c depends on direction, i.e. is a 
function on the unit sphere S' d_1 . We may then write V as r~ 2 (P(uj) + o(l)). 
Our first main result is 

Theorem 1. Consider the Schrodinger operator L = A+V in ~R d , where V is 
a bounded measurable potential equal to r~ 2 (P(uj) + t(r,uj)) for r > 1. Denote 
the eigenvalues of the operator A52 + P on L 2 (S d ~ v ) less than —\{d — 2) 2 by 
— /ii, . . . , — yU n . Assume one of the following two conditions: 

(i) The value —^(d — 2) 2 is not an eigenvalue of A S d-i + P on L 2 (S ,<i_1 ) 
and the function t(r,uj) is 0((logr) _( ^ 1+e - ) ) as r — > 00 for some e > 0, or 

(ii) t(r,uj) = 0((logr) _2_e ) as r 00, for some e > 0. 
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Then the counting function of L may be estimated as 



Nl(E) = vV* - (d - 2)74 + 0(1) (1) 

i=l 

In particular, if there are no —fa < —-Ad — 2) 2 then L has a finite discrete 
spectrum. 

Remark. The exponent in condition (ii) is sharp, in the following sense: there 
exists a potential V of the form of the form r~ 2 (— |(d — 2) 2 + i(r)) for r > 1, 
with t = 0((logr) -2 ), such that L = A+V^ has an infinite number of negative 
eigenvalues, and thus fails to satisfy ([IJ). See Section El for an example. 

In [HI], Kirsch and Simon proved (with the weaker error estimate 
o{\og{E~ 1 ))) in the special case that P is a constant function on the sphere. 

To prove our second theorem, we make additional assumptions. We work 
in three dimensions and assume that V is smooth, that for r > r > the 
potential V is exactly equal to P(Q)r~ 2 , and that the smallest eigenvalue 
— H\ is less than —1/4 while all others are strictly greater than —1/4. Under 
these assumptions we can prove 

Theorem 2. Let d = 3 and assume that the potential V satisfies the condi- 
tions of the previous paragraph. Let 

a = expf 2n V ( 2 ) 



Then there exists a G E stxc/i t/iat ; u>rf/i — £7 n t/ie nth eigenvalue of A + V 
counted with multiplicity, 



a as n — > oo. 



Moreover, if v n is the corresponding normalized eigenf unction, then v n is 
essentially supported in an annulus with radii proportional to a~ n ^ 2 , in the 
following sense: for any e > 0, there exist C + , C_ > so that for all n > 1, 
the L 2 mass of v n in the annulus 

C_a- n l 2 <r< C + a~ n / 2 

is at least 1 — e. 

The strategy of the proof of Theorem ^ is as follows. We introduce oper- 
ators L + and L_ which lie respectively above and below L, and then further 
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modify these to 'pseudo-Laplacians' L\ and L_ in order to eliminate diffi- 
culties near the origin. These pseudo-Laplacians L± have a domain which is 
a finite-dimensional perturbation of the domain of L±, in the sense that the 
intersection of the two domains is a closed subspace of finite codimension in 
each domain. It is a simple consequence of the min-max characterization of 
eigenvalues that the difference of the counting functions \N L i> (E) — N L± (E)\ 
is bounded by the maximum of these codimensions. Hence, to estimate 
Nl{E) up to 0(1) it suffices to get asymptotics for N L ^(E) and N L ^(E) 

which differ by 0(1) as E — > 0. The operators V± are arranged so that we 
can separate variables by introducing the eigendecomposition on the sphere, 
and we thereby reduce the problem to a very classical problem of estimating 
the number of zeros of the solution to an ordinary differential equation in 
the r variable. Our pseudo-Laplacians are somewhat analogous to operators 
used elsewhere in spectral theory to control difficulties caused by low-lying 
eigenvalues on a cross-section, as in [TT] . [S|, [2] for example. 

To prove Theorem we construct approximate eigenf unctions $„ for a 
sequence of approximate eigenvalues A„ asymptotic to ao~ n , where a is given 
by ©• For this to be effective we need the norm of the error term (H — A„,)$ n 
to be o(A n )|| ( I> ri ||2 as n —* oo. Our approximate eigenfunction <D n is equal to 
a zero- mode of A + V inside a ball of radius p n > tq, and an exact solution of 
the eigenfunction equation outside the ball, where we may use separation of 
variables, taking advantage of the assumption that V is exactly homogeneous 
in this region. In addition we need to add smoothing terms so that $ n is C 1 
across the interface. The approximate eigenvalues X n are determined by the 
criterion that the principal terms in the expansions of $ n inside and outside 
the ball -8(0, p n ) match in a C 1 way at the interface; this automatically makes 
the \ n asymptotic to a geometric progression (Lemma Hj). 

The result of the construction is that we can show that there is at least one 
eigenvalue in intervals of the form [aa n — o(a n ), aa n + o(cr n )]; note that these 
are non-overlapping intervals for large n. By Theorem Q however, there can 
be at most one eigenvalue in all but finitely many of these intervals, and at 
most finitely many eigenvalues lying outside the union of these intervals. We 
conclude that in some spectral interval [— Ao, 0) there is exactly one eigenvalue 
in each of these intervals and no other eigenvalues. It then follows that the 
$ n are close to true eigenfunctions, and the statement about the L 2 mass of 
the eigenfunctions can be verified by checking it for the $ n . 

In the final section we consider two examples. The first illustrates the 
remark after Theorem ^ The second example shows some limitations of a 
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heuristic approach of Fefferman and Phong to obtaining estimates for the 
counting function Nh(E) for an operator H . Their approach is to count the 
number of disjoint images of the unit cube in phase space W 1 x M n , under 
canonical transformations, which can be fitted into the part of phase space 
where the symbol of H is less than E. While this heuristic has been shown to 
give accurate estimates of counting functions (up to constants depending only 
dimension) in many cases [1], our example shows that no such heuristic 
involving connected regions of phase space will work in the present setting. 

2 Auxiliary Operators 

Our method is to estimate the operator L, and hence Ni(E), from above and 
below by two operators L± with a simplified potential. We then compare 
these operators with two further 'pseudo-Laplacians' L± which allow us to 
treat the 'near region' r < 1 and the 'far region' r > 1 separately. Let us 
choose a constant W such that | V(rr) | < W for all x G M d . Let 

t(r) =C(l + logr)-( 1+e \ 

where e is as in Theorem ^ and C is large enough so that \t(r, u)\ <t(r) for 
r > l 1 . We now define the potential operators V± by 

V±(x) = \x\- 2 P(uj) ±t(r), \x\ > 1; V±(x) = ±W, \x\ < 1 (3) 

and introduce the operators L± given by 

L± = A + V± 

with domain H 2 (R d ). Clearly we have 

L-<L< L + . (4) 

The pseudo-Laplacians are defined by the same formula as L± but 
with a different domain. The domain of V± consists of functions whose 
restrictions to the interior and exterior of the unit ball B d are in iy 2 2 (i? d ) 
and W / 2 2 (]R d /-B d ) respectively, together with some conditions at the boundary 
of the ball which depends on a choice of finite-dimensional subspace T of 
L 2 (S' d_1 ). Since the sphere S^ _1 of radius r centred at the origin is canonically 
identified with the sphere S^ 1 of radius 1, by radial projection, we can 
regard T as a finite dimensional subspace in L 2 (S^ -1 ) for any positive r. Let 

1 Under the assumption (ii), we will take t(r) — C'(l + logr)~( 2+e ) 
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n T , 11^ denote the orthogonal projections onto T, and onto its orthogonal 
complement, respectively, in L 2 (Sf~ l ). 

To describe these conditions, first note that ?p\ B d £ W^B 8 ) and ip\^d/ B d G 
W^iW 1 / B d ) implies, by the Sobolev embedding theorem, that the limits 
lim r jj ip(r, •), lim r ji d r tp(r, •), lim^i ip(r, •), lim^i d r ip(r, •) are all well-defined 
in L 2 (S d ~ l ). Indeed, i/j(r, •) is a C 1 ' 1 / 2 function of r with values in L 2 {S d ~ l ). 
The conditions we place on ip G Dom(L^) is that 

n T (lim^(r, •)) = n T (lim^(r, •)) = 0, 
n^(r, •)) G W 2 2 (R 3 ) 

Let us now prove 
Lemma 1. The operators L± are self-adjoint. 

Proof. From the definition of adjoint, we seek the operator with max- 

imal domain such that 

J(A + V±)M = j (6) 

for all ip G Dom(L±). If we restrict to ip supported inside B d , respectively 
M. d \ B d , we see that has a generalised Laplacian on these domains so 
<p\ B d G W 2 2 (S d ) and 0| Rd \ Bd G W 2 2 (M d \ We may then apply Stokes' 

formula to obtain 

/ / ^(A + V±)0 = lim(7 f^- / ^ 

V V ^Ti Vis? dn J S 2 dnj 

™ Vis* <9n J s? dn) ' 

As usual this means that the left and right hand sides separately vanish. Thus 
(£±)*0 = (A + V±)4>. To determine the domain, we consider the boundary 
terms at S^ 1 and write ip and as the sum of their projections into T and 
into T L . The cross terms cancel by orthogonality of T and T x , and we can 
analyze the terms in T and in T L separately. For the part in T -1 , the values 
of i/j and d n ip agree from r < 1 and r > 1, and this forces the values of <fi and 
d n ip to agree also. This implies that U^tp is in ^(M '). For the part in T, 
we have -0 = on both sides, but the values of d n ip are independent. This 
forces <p to behave likewise. This proves that Dom(L^)* = Dom(L±). □ 
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The domain if 2 (R 3 ) of L± can be described as those ip with iP\d 3 £ 
IVfl-D 3 ) and ^|k3/ D 3 e ^(IR 3 //) 3 ) satisfying the second condition in (jHJ, 
with the first replaced by 

n^lim^r, •)) = ^(lim^r, •)), 

U T (lim d r ?p(r, ■)) = nyflimcV^fr, •)) 

rj.1 rfl 

(in other words, the function Hr("0) and its normal derivative U(d r ip) are 
consistent across the unit sphere). This makes it clear that the intersection 
DomL-tflDomLj! is a closed subspace of codimension r = dimT inside both 
DomL-t and inside Domilf (with respect to the graph norm). It follows from 
this and from the minimax characterisation of eigenvalues that 

\N L± (E)-N L ^E)\<r. (8) 



Combining Q and (jHJ), we see that to prove (0) it suffices to prove that both 
L\ and Vt satisfy the asymptotic on the right hand side of (JTJ). 



3 Characterisation of Eigenfunctions of Ll 

We shall determine N L ^ (E) up to a constant by exactly describing all but a 

finite number of the eigenfunctions of L±. 

If ip is such an eigenfunction of U± with eigenvalue —A < 0, for r > 1 we 
may expand ip in the eigenfunctions {H i \°^ l of A S d-i + P(uj), 

oo 

4>(r,u) = '*rX i (r)H i (u), 

i=l 

and for r < 1 we may expand in the eigenfunctions of { Ji}^ of A S d-i ± W, 

oo 

^(r,u)=Y t Y i (r)J i (u). 
i=i 

The functions and {Ji}™ 1 are ordered by eigenvalue, and we recall 

the the eigenvalues of A S d-i + P{oS) are {—fMi}^. Let z/« denote the zth 
eigenvalue of A^d-i ± W. Because the potential V± is smooth in r for r / 1 
the Xi and Yj are smooth, and we may separate variables to obtain the 
following equations: 

r^X'l + (d- l)rX'i + (m =F t{r) - r 2 A)X; = 0, r > 1 (9) 
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r 2 Y" + (d- l)rY- - + r 2 \)Y,i = 0, r < 1. (10) 

Because ip e L 2 (IR 3 ) we must have Xi(r) — > as r — > oo. We claim that if T 
is chosen to contain enough of the functions and Jj, the eigenf unctions of 
Lj! are relatively simple. In particular, we wish to prove the following: 

Lemma 2. There is a k such that if T = span{Hi, . . . , Hk, J\ ■ ■ ■ , «/&}, aW 
but a finite number of eigenf unctions of V± are of the form 

k 

^{r,u;) = Y t X i (r)H i {u), r > 1, (11) 

i=i 

fc 

il){r } u) = ^2Y i {r)J i (u) } r < 1. (12) 

i=l 

Proof. We begin by recalling that the nonpositive spectrum of is discrete 
except for a possible accumulation point at 0. Because we may disregard as 
many eigenfunctions as we like, we may assume that the eigenvalue —A is 
restricted to an interval [— Aq,0], for some Ao > 0. Our strategy is to show 
that when Ao is sufficiently small, there exists k such that for % > k, the 
inequalities Fi(l)F/(l) > and Xi{l)X[{l) < are satisfied, with equality 
only if Xi or Yi are identically 0. 

This may be proven for any k with — > sup r>1 \t(r)\. If i > k, Xi(l) > 
and X'^l) > 0, and X[{r) is not positive for all subsequent r, there must 
have been some first r for which X-(r ) = 0. We know Xi(r ) > 0, and 
substituting r into (J0J) we obtain 

r 2 Xf(r ) = (-^ ± t(r) + r 2 A)X,(r ). 

We know — /ij ± t(r) > and r 2 A > 0, so this implies X"(ro) > which 
contradicts the assumption that X[ was becoming negative for the first time. 
Therefore X[, and so Xi, will be positive for all r, contradicting the as- 
sumption that Xi(r) — > 0. We may obtain a similar contradiction in the 
case Xi < 0, X[ < so we have Xi{\)X[{\) < 0. In the case X<(1) = 0, 
Xj(l) / we may show that X"(l) and X^l) must have the same sign, so 
Xi has a positive local maximum (negative local minimum) and we again get 
a contradiction. Therefore equality may only occur when Xj(l) = X^(l) = 0. 

We now wish to show that 1^(1)1^/(1) > for all but finitely many i. 
Because our eigenfunction is continuous at the origin, and each spherical 
harmonic (except for the constant function, which we may exclude by taking 
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k > 1) has average value 0, their coefficients in the expansion of if) must 
vanish at zero. 

Choose k so that u^, the kth eigenvalue of Agd-i ± W is greater than 
0. If Y%(1) > and Y/(l) < 0, then Y" must have a positive local maximum 
somewhere on the interval (0, 1). However substituting such a point into (fTU|) 
gives us Vi < —r 2 X < 0, a contradiction. The case Yi(l) < and Y({1) > 
is similar. Finally, if Y-(l) = then 3^(1) and Y"(l) will have the same 
signs, in which case we must again have either a positive local maximum or 
negative local minimum on (0,1). Therefore Yj(l)Y/(l) > with equality 
only if Yi(l) = Y({\) = 0, i.e. Y { vanishes on [0, 1]. 

To combine 1^(1)1^/(1) > and Xi(l)X^(l) < and obtain a contradic- 
tion, we observe that 

/ n^M-V (!3) 
Jan dn 

is the same when the derivative is taken from the inside and outside of the 
unit sphere. This follows from the fact that U^ip e (R d ). Expanding in 
eigenf unctions from the inside and from the outside then gives 

oo oo 

o > x<0)x'iQ) = E WWO) * 

i=k+l i=k+l 

Therefore X;(1)X;(1) = and Y;(1)Y/(1) = for all z > k, so only the 
spherical functions contained in T may be present in any of our remaining 
eigenf unctions, proving the lemma. □ 



4 A Sturm-Liouville Reformulation 

Lemma 1 effectively reduces the study of the discrete spectrum of Li± to the 
finite number of one dimensional Sturm-Liouville problems (|11|) . (J12)) . We 
may rewrite Q and (JTUJ) as 

r w x, = xv + d -^x[ + ^±^x = \Xi (14) 

and 

n^Yi = Y" + — Y/ - V -\Y t = XY. (15) 

The Xj have boundary conditions Xj(l) = (due to the domain condi- 
tion (0)) and limXi(r) = 0. Yi must satisfy Yi(0) = Y(l) = for % > 1 and 

r— >oo 

Yi(l) = 17(0) = 0, the Neumann condition at again arising from smooth- 
ness of if) at the origin. To analyse the spectrum of T and Q we use the 
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following well-known link between the number of eigenvalues of the operator 
less than E and the number of zeroes of a solution to (|14|) with A = E; 
see for example Theorem XI 1 1. 8. 

Theorem 3. If n(T^, E) is the number of zeros of the solutions ofT^X = 
EX, X(l) = for E > 0, and N(T^, E) the counting function of the operator 
T M with boundary conditions as above, then N{T E) = u(r w , E) for E > 0. 
The analogous result holds for the Q Ui . 

As we are considering only a finite number of operators, the 0(1) may be 
taken to be uniform. It remains to estimate the number of zeros of a solution 
to T^X = EX, X(l) = 0. To do so we transform by substitution: setting 
h(r) = Xr-( d ~ 2 ^ 2 we obtain the equation 

r 2 h" + rh! + {fi- -{d - 2) 2 - r 2 E + t(r))h = 0, 

and setting g(r) = h(e r ), 

g" + (/i - i(d - 2) 2 - Ee 2r + t(e r ))g = 0. (16) 

Proof under hypothesis (i). In this case, by assumption there is no eigen- 
value equal to —\{d — 2) 2 so we can split our analysis of (fTB^) into the cases 
— /i < — -Ad — 2) 2 and —ji > —\{d — 2) 2 . In the first case, normalise the 
coefficient a — \{d — 2) 2 to be 1 by letting a = , 1 = and defining 

r 4V ' v / Ai -(d-2) 2 /4 

/(r) = g (ax). The final form of our equation is then 

/" + (l - a 2 Ee 2ar + t(e ar )) / = (17) 

We may think of this equation as describing paths (/, /') in a phase plane, 
so that estimating the number of roots becomes a question of estimating the 
rate at which our point travels around the origin. Let 9 be the angle between 
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the point (/, /') and the positive x axis. We have 

f -if 
= tan — 

1 f'f - {f? 

i + (W) 2 ' P 

f'f - if) 2 
p + if) 2 

(1 - a 2 Ee 2ar + t(e ar ))f 2 + (f'f 

P + (f? 
(a 2 Ee 2ar - t{e ar ))f 2 

' ' P ■ iff 

< a 2 Ee 2ar + t(e ar ) 

Therefore for small r, the point (/, /') will travel clockwise around the origin 
at a nearly constant rate. The error term a 2 Ee 2ar +t(e ar ) becomes very great 
for large r, but this is not a problem because as soon as 1 — a 2 Ee 2ar + t(e ar ) 
becomes negative the character of the ODE (fTTj) changes. In particular, if tq 
is such that 1 — a 2 Ee 2ar + t(e ar ) < for r > ro, then by (fTTj) / cannot have 
a positive local maximum (or negative local minimum) with r > ro and so 
can only have one root greater than ro- Therefore to estimate the number of 
times (/,/') winds around the origin we only need to consider the behaviour 
on the interval [l,ro] where 1 + ||t||oo > a 2 Ee 2ar holds. For all E 

f a 2 Ee 2ar dr < -aEe 2ar ° < 1 + l|t||o ° , (18) 
Jo 2 2a 

so the contribution to the error from the a 2 Ee 2ar term is bounded. The form 
of t(r) implies that J °° t(e ar )dr < oo giving us 

0(r o )-0(O) = -r„ + 0(l), 

so the number of roots of / is r /7r + 0(1). Calculating r in terms of E 
(up to a 0(1) term) gives r = log(l/E)/2a + 0(1). Therefore we have the 
number of roots equal to j-y/fi — l/4hi(l/E) + 0(1). (We also remark that, 
since there are only finitely many eigenvalues below —d(d — 2) 2 /4, there are 
only a finite number of values of a and so the 0(1) estimates are uniform in 
a.) 

In the case — 1/4 < — /x, we have 1/4 — /x > and so may rewrite our 
equation for g as 

g "-(l/A-fi + Ee 2r -t(e r ))g = 0. 



e 

de 

dr 
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Because t = o(l), there will be an r independent of E such that (1/4 — ji + 
i?e 2r — t(e r )) > for r > r and E small, and after which g can have at 
most one more root as discussed above. Also, the number of roots that g 
may have in the interval r < r$ is bounded for small E by the continuity of 
solutions of ODEs in their coefficients in compact domains. Therefore when 
—1/4 < —/ji, the number of roots is bounded as E —* 0. 

A similar analysis of the operators f2j is much easier, because they act on 
functions with compact domains and so solutions to QiY = EY, Y(l) = 
are uniformly continuous in E. Therefore the number of roots of solutions is 
bounded as E — > 0, so these problems only have a finite number of positive 
eigenvalues. 

We therefore have 

w) = ^i;v^i+o(i) 

Z7T * — ' 

i=l 

as required. The reverse estimate on the potential, and the counting function, 
is performed in exactly the same way. 

Proof under hypothesis (ii) In this case, we have to analyze equation (jTSjl 
also in the case that — /i = —1/4, but now we take that t(r) = C(\ogr)~( 2+e \ 
Let T(r) = t(e r ), so that \T(r)\ < Cr-( 2+€ \ We have to show that the 
number of zeroes of the solution of (JTfij) is 0(1) as E — > 0. 

Since the —Ee 2r term is negative, by the Sturm comparison theorem it 
suffices to show that the number of zeroes of the solution of 

g" + Tg = (19) 

is finite. Also, by the Sturm comparison theorem, we may assume that T > 
(otherwise, we can consider max(T, 0)). 

To do this, let q be a root of g, and assume (by scaling g) that g'(q) = 1. 
(We cannot have g'(q) = unless is the trivial solution.) To show that g has 
no root larger than q, it suffices to show that 

POO 

/ -g"(r) dr < 1, 

Jq 

since that will show that g is increasing on [q, oo). We have 

roo roo roo f < 

/ -g"(r)dr< Tgdr<C r~ {2+e \r - q) dr < -q~\ 

J q J q J q ^ 

This is indeed less than 1 for q sufficiently large. Thus for ro sufficiently 
large, depending on C and e, we see that g can have at most one root on 
[ro, oo). The number of roots in the interval [l,r ] is uniformly bounded 
for E G [0, Eq] so this gives a uniform bound on the whole real line. This 
completes the proof of the theorem under hypothesis (ii). 
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5 Preliminaries to the proof of Theorem [21 

Our method of proof is the technique of approximate eigenfunctions, de- 
scribed in the following well-known lemma: 

Lemma 3. Let A be a self-adjoint operator on a Hilbert space H . Suppose 
there exists u ^ G H with 



Then there exists a point of spec(A) in the interval [A — e, A + e] . 

Now let /i > e, and suppose that the interval [A — //, A + //] contains only 
one point of spec(A) which is an eigenvalue of multiplicity 1. Let v be a 
normalized eigenf unction for this eigenvalue. Assume also that (u,v) > (if 
not, then replace v by av where a & C, \a\ — 1 is chosen so that (u, av) > 0). 
Then if u satisfies (J2U|) . 



Proof. If there is no spectrum of A in the interval [A — e, A + e] then, since 
the spectrum is closed, there is e' > e such that there is no spectrum of 
A in the interval [A — e', A + e']. By the spectral mapping theorem, then, 
\\(A - A)" 1 ! < 1/e'. This contradicts flU. Hence spec(A) n [A - e, A + e] 
must be nonempty. 

To prove the second part of the lemma, we write u = cv + w where 
c = (u, v) > and w is orthogonal to v. Then 



\\(A- X)u\\ < e\\u 



(20) 




\\(A-X)u\\ > \\(A-X)w 



\\(A — X)cv\\ > fi\\w\\ — ce 



which implies that 




2e 



Now 



u — v\\ < \\w\\ + \\v 



cull < h 1 — C. 



Finally 




so 



u — f || < 



6e 



□ 



13 



We therefore wish to find a sequence of approximate eigenvalues — A n , 
tending to zero, and functions $ = $ n such that || (L + A n )$„,||/||$ n || is small. 
Because V is unperturbed for r > tq, we may separate variables to construct 
exact solutions of L$ = A n $ in this region. However, these exact solutions 
are very poor approximations where V is perturbed. One way around this 
would be to use a given eigenfunction of of L in this region, however if we 
choose an eigenfunction \l/ with eigenvalue — v < for this purpose we have 
(L + A n )$ = (A n — in the region where \1/ is used. This gives us an error 
with a | A n — u\ term in it which remains rougly constant as A n — > 0, which 
is bad from the point of view of applying Lemma El A better choice for \1/ is 
a well behaved zero-mode of L, because now (L + A n )$ = A n \l/ and we have 
an error behaving as A n . The existence of such a function is guaranteed by 
the following result: 

Proposition 1. Let the operator A + V in IR 3 be as described above. Let 
— > 1/4 be the smallest eigenvalue, and J\{oS) the corresponding eigen- 
function, of the spherical operator A52 + P. Denote by N e denote the null 
space of A + V acting on (x)~~ e L 2 (R 3 ), where (x) = (1 + |x| 2 ) 1//2 . Then 
for small positive e the dimension of iV_]_ e /iV_i4. e is one dimensional. This 
quotient space is spanned by a function ^ satisfying 

(A + V)ty = 

and having the asymptotic expansion 

^ = r - l/2 cos(vVi - l/41n(r) + c) Ji(w) + 0{r~ 1/2 ~ 5 ), r -> 00 
for some c G R and some 5 > 0. 

Proof. These results are implied by the relative index theorem in chapter 6 
of jT2] • Let us write r for a function that is > 1 everywhere on M 3 and equal 
to I a; I for \x\ > po- We first note that we can write 

A + V = r~ 5 / 2 ( - (rd r ) 2 + A S 2 + P + r 1 ' 2 

for r > po- Therefore the null space of A + V on (x)~ 1±e L 2 (]R 3 ) is equal to 
the null space of an operator Q on r ±e L 2 (lR + x S 2 ;fib) where Q takes the 
form — (rd r ) 2 + A52 + P + | near r = 00 and /j.^ is a smooth measure taking 
the form doodr/r for r > p (where du is the standard measure on S 2 ). The 
operator Q is an elliptic b-operator in the sense of ^2] • By the relative index 
theorem of Q is Fredholm as a map 2 

Q : r- t H 2 (R + x S 2 ; p b ) -> r- e L 2 (M + x S 2 ; /i 6 ) (21) 

2 Here i? 2 is the b-Sobolev space of order two, given by the functions whose derivatives 
with respect to angular derivatives <9 W and with respect to rd r up to second order are 
square integrable. 
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for all e 7^ in a neighbourhood of 0, and, with Nf denoting the null space 
and ind e Q the index of Q acting as in ([21)1 . for small e > 0, 

ind_ e Q — ind e Q = 2. 

(This number is equal to the number of indicial roots, i.e. complex numbers 
r for which A52 + P + 1/4 + r 2 is not invertible, with imaginary part in the 
interval [— e, e]. Here we have r = ± a//Ui — 1/4 as the only such indicial roots, 
provided e is sufficiently small.) Since Q is self-adjoint on L 2 (R + x S 2 ;/i b ) 
and by elliptic regularity we have 

ind e Q = dim Nf - dim N Q £ . 

Therefore ind e Q = — ind_ e Q, and we have for sufficiently small e > 

ind_ e Q — 1. 

This implies that dim(A^„ 1 _ e /A^_ 1+e ) = 1. If we take any function \l/ in 
N^i_ e \ N_i +e , then the regularity results of [12], e.g. Proposition 5.21, show 
that \l/ has a conormal asymptotic expansion as r — > 00, i.e. an expansion in 
powers and powers of logarithms of r with coefficients that are eigenfunctions 
of the operator on the boundary. This implies the asymptotic property stated 
in the theorem. □ 

Our approximate eigenfunctions, constructed in the following section, will 
be built out of the function \l/ from Proposition ^ close to the origin, and 
an exact solution to the eigenfunction equation in the far region, where we 
can use separation of variables. If the appropriate spherical eigenfunctions 
are {J i \°l l with eigenvalues {—f^i}^, we may expand F as 

00 

F(r,u) = Y t f i (r)J i (u), 

i=l 

where for r > the components fi obey the differential equation 

r 2 /r + 2r/;+(/i i -r 2 A)/ l = 0. (22) 
Letting fi = r~ 1//2 /j, the function fi satisfies the modified Bessel equation 

r 2 f:' + rf: + (^-\-r 2 X)f t = 0. (23) 

There is a one-dimensional space of solutions which are exponentially de- 
creasing as r — > 00, given by multiples of the MacDonald functions fi = 
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K r^z^ (V^ r ) ( see section 9.6). It is a standard fact that these func- 
tions K v {r) have no zeroes for r > (pQ, p377). 
Let us define the exponents a«, /3j by 

a, = ^ - f^ . ft ^ 1 + f^ . (24, 
Then the functions r°% r* 3 * solve (J22|l with A = 0. Let 

/oo 
(p 2 + l) Q *e- ipr dp. (25) 
-oo 

Then Xi is a multiple of r" x l 2 K ^/ 1 y 4 _ M ( r ) an d solves (|22p. For each % > 2, 
let 

X, A (r)=X,(rv / A); (26) 

this is an exponentially decaying solution of (|2*2*|) for arbitrary A, with the 
explicit expression 

/oo 
(p 2 + \) a >e- ipr dp. (27) 
-oo 

We now obtain some estimates on the needed in the proof of Theo- 
rem |21 The identity (J2*7j) lets us estimate the value of Xf for < —1/2, i.e. 
i > 2, by 

POO 

X*(r) = (rV\) a > / (p 2 + l) ai dp (l + o a .(l)), n/A->0. (28) 

J — oo 

Here the o term is uniform as % — > oo, but blows up as — > —1/2 (which 
corresponds to /ij — > —1/4); however, since there are only a finite number 
of OLi in this range, for a given potential, we may take this estimate to be 
uniform in i > 2. This implies that 

Xi(r) > C(rVX) ai \2ai + l]" 1 for sufficiently small rVX. (29) 

For CKj < —3/2 we may use (|27p to prove 

^^ = ^(l + 0((rv^) 2 )), (30) 

where the O term is uniform as i — ► oo. 

On the other hand, series solutions at r = show that if —1/2 < ol% < 
—3/2 and a« is real, 

^^ = ^(l + Q! (ln(rv / A)(rv / A)- 1 - 2 -)), rV\ - 0. (31) 
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Here the O ai is not uniform as ctj — > — 1/2, but as above, since there are only 
a finite number of a, in the range [—3/2, —1/2) for a given potential, we may 
take this estimate to be uniform in i. If we define 



K = min (2, min( aA) > 

v ' i>2 v 2 

then we may combine the two estimates (JBTij) and (j3*Tj) in the form 



(32) 



X A fr) 



— (l + 0((n/A)*)), «>2, rv/A^O. 



In the case i — 1, there is a constant d such that 
X A (r) = r - 1/2 A" 1/4 (cos(V)"i-l/41n(rVA) + d) + O(rVX) 



(33) 



(34) 



and 



(X*)'(r) = r- 3 / 2 A- 1/4 ( - l/2cos( v / /ii - l/lln{rV\) + d) 

- 1/4 sin( V/zi - 1/4 ln(rVA) + d) + O(rVX)) , (35) 



where both functions 0(ry/\) are smooth away from 0. 

We also need estimates on L 2 norms. Using ()27|) we may bound the L 2 
norm of X A Jj restricted to M 3 \B(0,p) as follows: 



I Y x T II 2 

J i\\L 2 (R3\B(0, P )) 



(X i (rv / A))V 2 rfr 

/oo / /"oo \ 2 

(rV\) 2a * \J {p 2 + l) a 'e- iprVX dpj r 2 dr 

/oo / poo \ 2 

(p 2 + l) a< e- ipn/X dpj d(rVX) (36) 





-1/2 


< p 2Qi+2 A ai 


-1/2 




-1/2 



{p 2 + l) a >e~ ipr dp 



\Oi\\2 



< o ai (i)p 2at+2 x at - 1/2 



Here the 0(1) is uniform as ctj — > — oo, but not as ai — > —1/2; again, as 
above, for a given potential, we may take this estimate to be uniform in i. 
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6 Approximate eigenfunctions 



Our approximate eigenfunctions, for suitable values of the eigenvalue —A, 
A > 0, will be built out of three components. These are the zero mode \I/ 
from Proposition d close to the origin; an exact solution to the eigenfunction 
equation far from the origin; and smoothing terms supported in an interme- 
diate region, which ensures that the two pieces fit together so as to lie in the 
domain of L. To define this, let h(t) be a function which is supported on 
the interval [1, 2], equal to t — 1 on [1, 1.5], and smooth on (1, 2]. Hence h is 
smooth except at t — 1 where it has a jump in the first derivative. For the 
nth approximate eigenvalue, — A n , we will choose a radius p n , growing as a 
negative power of A„, and define 

r<p n , 

$ ~ " ~ [EZi (&Xt(r) + xMr/Pn))jM r > p n . (3?) 

Here fa is a coefficient chosen to ensure the continuity of $ at r = p n , and 
the coefficient Xi is chosen so that the radial derivative of $ is continuous at 

r = Pn- 

We may perform separation of variables on \1/ in the region r > tq to 
obtain an expansion 



^(r,u) = Y i Y i (r)J i (u). 



8=1 



The functions Yi must satisfy ()22|) for A = and r > tq, which has general 
solution Yi(r) = ar ai + br^ where a^,/?, are defined in (|24j) . For % > 2, on 
and Pi are real and aij < —1/2 < (3^ since ^ = 0(r -1 / 2 ), we have 

Y i {r) = ^j i r a \ r>r (38) 

for some coefficients ipi. 

We choose A n so that the smoothing term for i = 1 is not needed, i.e. so 
that we can choose Xi — 0. We therefore need to choose A so that 



(X^(p)Y 1 (p) = Y;(p)X^ P ). (39) 



Lemma 4. Let 



, -2nn — C . 
exp ( , | (40) 

^i - ! 

for some CeR and let p n be a sequence of real numbers > r with p n y/C^ — > 
0. Then there exists a constant C and a sequence (A n ) ; satisfying J3[PJ) at 



18 



p = p n which are asymptotic to £ n in the sense that 

= 1 + 0(p n \fQ) as n -> oo. (41) 

Remark. All the O(-) and o(-) estimates in this section are uniform as i — > oo, 
A — ► and p — ► oo provided p\/A — > 0. 

Proof. Since we know that 

Yi(r) = r -1 ^ 2 cos(a/ p\ — 1/4 log r + c), r > r , 

we have 
Y[{r) 



_ 3 
/■ 2 



Q cos(a//Ui - 1/4 log r + c) + vVi ~ V^si^y^ - 1/4 log r + 



r 2 cos( a/ /ii — 1/4 log r + c — 0') 



where cos#' = (2 V //I7) Similarly we have 

X${r) = r- 1/2 A- 1/4 (cos(VA*i - l/A\og{rV\) + d) + 0(rV\) 
and 

j-jtf(r) = -r- 3/2 A _1/ Vi(cos(v^ - l/41og(rv / A) + d - 9') + 0(rV\)^ 
as ry/X 0. We want to solve 

Y l {r)±-X^r) = Y{{r)X^ r ) 

which is equivalent to 
cos( \JpT 1 - 1/4 log r + c) ( cos(\/p 1 - 1 /4 log(r v^X) + d - 0') +0(tVA)J 

= cos( v / /Ui - 1/4 log r + c-9') (cos^^i - l/41og(rV^) + d) + 0(tVA) J 

at r = p n . Using trigonometric identities this reduces to 

sin & sin {\J p\ - l/41og\/A + d - c) = 0(p n VX). 

The sin^' can be neglected here since this is a positive quantity depending 
only on p,\. We note that the 0(p n V~X) term on the right hand side is a 
continuous function of A. Therefore, by the intermediate value theorem we 
get a sequence of solutions of the form 

1 , , c — d — irn . rr— 
-logA n = +0(p n \/\ n ) 

* y/Hl - 1/4 

for n sufficiently large, which satisfies (}4*T]) with C = 2(d — c). □ 
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We now analyse the matching conditions required at 3 r = p. In order to 
have continuity of $ at r = p we choose, for i > 2, 



(42) 



and in order to have continuity of <9 r $ at r = p we choose 

Xi = p(Yl{p)-HX})'{p)). 

(We remark that Xi has no zeroes for i > 2, so the denominator in (}4*2*|) is 
never zero.) Using ()42j) . (jHHj) and (}3*3*j) we find that 



Xi = Wi^iP - fap A 

Xi\P) 

Ap ai (oa - P{ ^ P) ) = ^0(\a t \(pV\) 



(43) 



where K is defined in 

Lemma 5. The expression (13 7|) converges in the graph norm on the domain 
ofL. 

Proof. Since \l/ is smooth, it is enough to show that $ is in the Sobolev space 
VFf (M 2 \ 5(0, r )). In this set we have $ = ^ Sj, where Sj is given by 

il)iYi(r)Ji(uj), r < p 

iX t x {r) + xMr/p))Ji{uj), r > r . 

Due to our matching conditions, each Si has two derivatives in L°°, and Si is 
rapidly decreasing together with all derivatives as r — > oo. Therefore, each 
term Sj is in the domain of L. We have to show that the sum converges in 
the graph norm. Since \1/ G W2{B(0, p)) and = 0, it is sufficient to show 
convergence of the infinite sums 

^2 / (0.X* + Xl h P ) 2 r 2 dr (44) 

i=2 J P 

and 

r 2p 



r p 2 /// \ 

£ || (L + A)(<^ A + x.M^f = E / + + B - A ) 

i=2 i=l r 

(45) 



5 For ease of notation we will drop the subscript n from X n and p n from now on. 
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where we write h p (r) = h(r/p). We may further reduce (J44j) to showing the 
convergence of 

oo /»oo 00 /*oo 

E / ^) 2 r 2 dr and E / X \h 2 p r 2 dr. (46) 

i=2 i=2 J P 

Similarly, ()45|) reduces to showing convergence of 

" 2P ' { tun , fy^\2 2 > /" 2P (XiP>ihp\2 2 



r^P h' rip t. 

£/ (»(*;+ 5?)) v*. W (^v* 

i=2 J ^ i=2 J /> 

and E] / (\xihp) 2 r 2 dr 



(47) 



(recalling that /i p is supported in [p, 2p]). We may obtain useful bounds on 
the sums in terms of H^mll where ^2r is the restriction of \1/ to the disk of 
radius 2r , using the estimate 

p2r 

/ r 2a *r 2 dr > C2 Ql r 2Ql (48) 

Jro 

for some constant C independent of i. Equation ()48|) gives us 

00 r2r 



00 > n^2r ir > E / ^ 2r2Qir2 ^ > c E 2ai ^ r o ai 

i=l j=l 



(49) 



This gives us in particular the convergence of the right hand side, and this 
allows us to show the convergence of the second component of ()44j) . 
If n is such that a n < —3/2 then we can then estimate, using 



00 /»oo 00 /»oo 

E / (X*h p ) 2 r 2 dr< J2^P 2a ^O((pV\) 2K ) h 2 p r 2 dr 

-i+l J P i=n+l J P 

00 

= E ^^((pv^) 2 *) 

i=n+l 

00 

< E ^(2r ) 2 - +3 a 2 O((pv / A) 2 ^), (50) 



i=n+l 

00 

<0((pv / A) 2 ^)sup(2^a 2 ) E V^r** 

i=n+l 

< 0((pv / A) 2 ^)sup(2^a 2 )||vI/ 2ro || 2 by © 

= o((pVa) 2 *) 
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which shows convergence of the second term in (|46|) . 

A similar estimate proves convergence for the first and third terms in 
(14 7|) . As for the second term, we need to take into account that — > oo. 
Using supj 2 ai afpf < oo (which follows from Weyl asymptotics for the p^, 
we can apply the same argument to show convergence of the second term of 

dm). 

It remains to show the convergence of the first term of (J46|) . We have 



OO pQQ OO 

£ / {faXtY^dr < 0(1) 0?P 2 " l+2 A a ^ 1/2 by © 

=n+l P i=n+l 

= ow E dS)) 2 " 2 ' > ' +2A ' >, " /2 by 



0(1) (pv/A)^ 20 ^^ 4014 " 2 ^ ^ 172 by (J2 



(51) 



i=n+l 
oo 



= 0(1) £ ^a- 1 / 2 ^ 2 

■t=n+l 

= 0(A-^ 2 ) 

using the same reasoning as in (|5L))1. □ 

In order to apply Lemma El we need accurate bounds on the norms of $ 
and (L + A)$. 

Lemma 6. Assume that p = \~ l / 2+5 / 2 , where S > is sufficiently small. Let 
$' be $ — 0iX^Jx restricted to M? \ B(0, p). Then for some e > (depending 
on 5) we have estimates 

l^lh^cA- 1 / 2 , (52) 
||$'|| 2 = 0(A~ 1/2+e ) (53) 

and 

||(L + A)$|| 2 = 0(A 1 / 2+e ). (54) 
Proof. To obtain a lower bound on ||$||2 it is enough to estimate 

||0iXl Ji\\l 2 (m 3 \b(o,p) 
in view of orthogonality of the Jj on the sphere. By scaling, we compute 



{X^{r)fr 2 dr = A~ 3/2 / A 2 (r)r 2 dr 

JpVx (55) 

> cA~ 3/2 since pVX 0. 
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On the other hand, we may estimate (j>i by combining the equations 



( p 1 X?(p)=Y 1 (p), 
ct> l {X x l )\ P )=Y{{p) 

to get 

ti{X X M 2 + P 2 (*i A )'(p) 2 ) = Y 1 (p) 2 + p 2 Y[{pf. 

From (f3~4"j) and (f33J) we see that the left hand side is bounded above and 
below by nonzero constant multiples of 2 p _1 A~ 1//2 and the right hand side is 
similarly estimated by p -1 (note that ip\ 7^ so the lower bound is nonzero), 
so 0i ~ A 1//4 . Therefore the norm of <piX^Ji behaves as A _1//2 . 

To prove we note that (joHj) and (joTj) show that if we take the sum of 
Si over i > n + 1 instead of z > 2 then we get a bound 0(A~ 1//2 ) for the square 
of the L 2 norm. To bound the terms with 2 < i < n we can adapt these 
estimates. In the second line of (p)Uj) . if we remove a factor of p 2 then the 
rest of the argument follows, since we have 2a, + 1 < for i > 2. Therefore 
we can end up with 0(p 2 ) for the square of the norm of each s,, which is 
0(\~ 1+s ). In (J51|) we can reach the second last line, and then we observe 
that for 2 < i < n, the factor p 2oti+2 is bounded by p 2a2+2 , giving an overall 
estimate of 0(A~ 1//2 p 2a2+2 ) for the square of the norm of Sj. This is 0(A~ 1+2e ) 
for e = (— 2«2 — l)/4 > which is also of the claimed form. 

To prove (p>4"|) we must estimate the sums in (@7J). To do this we use the 
computation ()50|) . modifying it as needed. For the first and second terms in 
(|47|) we have an extra factor of p -4 , as compared to the computation in (|5U|). 
so the sum of these terms (for i > n + 1) is 0(p~ 4 ). In the case of the third 
term we have an extra factor of A 2 , so this term contributes 0(A 2 ). 

Of course we have to bound the terms where i < n. Since there are only 
a finite number of these terms we can ignore the pi and ip^. Thus, for the 
first and second type of term in (J47|) . for a single i < n, we get a bound (from 
the second line of (J5U)) . recalling we have an extra p~ 4 ) 0(p 2ai ~ 1 (p\/\) 2K ) = 
0(p 2a2 '\pV\) 2K ). For the third term we get 0(A 2 p 2a2+3 (pv / A) 2X ). In sum- 
mary, we obtain 

|| (L + A)$|| 2 = 0(p- 4 ) + 0(A 2 ) + 0(p 2a ^ l {p^X) 2K ) + 0(X 2 p 2a2+3 (pV\) 2K ). 

Since p\^X — > and K > this is 0(p~ 4 ) + 0(p 2 ° 12 ^ 1 ). Now we choose 
p\/A = a/A , i.e. p = A^ -1 / 2 )/ 2 where < 5 < 1/2 is chosen so small that 

p- 4 = 0(A 1+2e ) and p 2 ^- 1 = 0(A 1+2e ). 

This requires 5 < min(l/4, — 02(1 — 2a2)~ 1 ). With this choice of p as a 
function of A we have ||(L + A)$|| 2 = 0(A 1+2e ) and hence (JSj). □ 
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7 Proof of Theorem [21 



Applying Lemma|3]to $ An , we see that estimates (J52)l and ()54|) imply there is 
a point of spec(L) in the interval [— A n — 0(A* +e ), — A n + 0(A* +€ )]. Combined 
with Lemma EJ this shows that for all n > no 

spec L n [-£„ - 0(^ +e ), -£„ + is nonempty (56) 

where £ n and no are as in the lemma. As discussed in the introduction, 
we may combine this with Theorem^ to show that in some spectral interval 
[— Ao, 0) there is exactly one eigenvalue in each of these intervals and no other 
eigenvalues. Consequently if — E n is the n th eigenvalue of L, counted with 
multiplicity, then 

,. E n 
hm — 

converges. This is equivalent to saying that if 

2tt 



a = exp 



y/JTi - 1/4 

then the limit 



hm ^ (57) 



n- 



converges, which proves the first part of Theorem |21 

We now wish to apply the second part of LemmaElto show that if v n is the 
normalised eigenfunction corresponding to —E n e [— Ao,0) then ||$ A ™ — v n \\ 
is small. Here 4> A " = $ A "/||$ An ||2 is the L 2 - normalized function. We have 
just seen that for n and C sufficiently large there is exactly one eigenvalue 
in the interval [— £ n — C£ 1+e , — £ n + CC, 1+e ] and no other eigenvalues in the 
interval [— Ao,0). Equivalently, there are integers no, n\ such that 

E n G [£ n+ni - Cf^+^fn+m + CCn+nil for a11 n>n . 
Lemma El gives 

||<l A ' i+ ' l i - v n \\ = 0(d) for all n > n . (58) 

By (f52J) and (JH3J), the L 2 mass of $ An is essentially supported outside B(0, p n ) 
by the first spherical eigenfunction Jj. Moreover, since |^| < cr -1 / 2 , the L 2 
mass of <3> An in B(0,p n ) is 0(p n y/X^) ~~ > 0. Finally, because X A is a scaled 
version of the fixed function X\, and since X n /cr n converges by (|4ip. there 
exist C + and C_ such that 1 — e of the mass of X An J t is supported in the 
annulus 

C_ a - n/2 <r< C+a- n/2 . 



24 



We conclude that there are C + and C_ such that 1 — e of the mass of $ An +"i 
is supported in the annulus 

C _ a -n/2 < r < C +(J -n/2 

Combined with ()58|1. this proves the second part of the theorem. 

8 Examples 

We now give an example illustrating the remark following Theorem ^ That 
is, we show that there is a t(r) with t(e r ) = 0(r~ 2 ), such that the Schrodinger 
operator with potential equal to r~ 2 (— \{d — 2) 2 + t(r)) has an infinite number 
of eigenvalues. 

Following the reasoning of the previous few paragraphs, it is enough to 
exhibit a function T{r) which is 0(r~ 2 ) such that the ODE (|19|) has an 
infinite number of zeros on [1, oo); then the number of zeros of must 
grow without bound as E — > since, on each compact set, the solution will 
converge uniformly as E — > 0. 

We choose a C 2 function g(t) on the interval [1,4] so that g(l) = 1, 
g(2) = — 1 and g(4) = 1, such that g is strictly monotone on the intervals 
(1,2) and (2,4), and so that g" vanishes on a neighbourhood of the point 
in (1,2) and the point in (2,4) where g(t) = 0. Finally, we assume that 
g'(l) = 4g'(4) and g"(l) = 16g"(4). It is clear that such a function exists. 
We now define g on [1, oo) by requiring it to be multiplicatively periodic with 
period ratio 4, i.e., that g(4t) = g{t). The function g(t) is now C 2 on [1, oo). 
The ratio 

Tit) = 

is also a C 2 function since we required g" to vanish near where g = 0. By 
the periodicity assumption, T(t) = 0(t~ 2 ). Clearly g solves the ODE (fT^j) 
and has an infinite number of zeroes. 



Our next example is related to the work of Fefferman-Phong. We work in 
dimension 3 for convenience, though the construction works in all dimensions 
d > 3. Let (9, 4>) be the usual coordinates on S* 2 , in which the upper hemi- 
sphere is {9 < 7r/2}. We consider a potential function P{9) on the upper 
hemisphere of S 2 which is equal to —1/3 for 9 < n/2 — 2e, equal to zero 
for tt/2 — e < 9 < tt/2 and monotone in between. Let P ev denote the even 
continuation of P to the sphere, and P Q dd the odd continuation. We assume 
that e is sufficiently small (less than 0.01, say). 
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Lemma 7. (i) The lowest eigenvalue of the operator A + P cv on L 2 (S 2 ) is 
less than —1/4. 

(ii) The lowest eigenvalue of the operator A + P G dd on L 2 (S 2 ) is greater 
than -1/4. 

Proof. To prove (i), we note that the quadratic form 

Q cv (f)= [ \vf\ 2 + p cv \f\ 2 dio 

Js 2 

takes a value — 1/3 + 0(e) on the L 2 -normalized constant function, which is 
< — 1/4 for e sufficiently small. 

To prove (ii), we take an arbitrary / e L 2 (S 2 ) with ||/|| 2 = 1, and write 
it / = f ey + /odd in terms of its even and odd parts with respect to reflection 
in the equator 9 = tt/2. Then, due to the symmetries of the sphere, 

/ \Vf\ 2 du= [ \Vf cv \ 2 du+ [ \Vf odd \ 2 duj. 
Js 2 Js 2 Js 2 

On the other hand, / Q dd is in the domain of the Dirichlet Laplacian on the 
upper hemisphere. The lowest eigenvalue of the Dirichlet Laplacian on the 
upper hemisphere is 2 (this follows since eigenfunctions of the Dirichlet Lapla- 
cian on the upper hemisphere are the restrictions of odd spherical harmonics, 
and the smallest eigenvalue of an odd spherical harmonic is 2). Hence 

/ |V/ odd | 2 ^>2 / \f odd \ 2 duj 
Js 2 Js 2 

and so 

/ |V/| 2 ^>2 / \f odd \ 2 doo. 
Js 2 Js 2 

Thus 

Qodd(/)= / |V/| 2 + P od d|./f^ 
Js 2 

>2||/ odd ||2+ / -P dd(|,/ev| 2 + I /odd 1 2 + 2Re/ cv / dd) du. 
Js 2 

Since P odd is odd and |/ cv | 2 , | /odd 1 2 are even, these terms vanish and we are 
left with 

W) >2||/ odd || 2 -2 / |Podd||/cv||/odd|^ 

Js 2 

> 2 1 1 /odd 1 1 2 _ g 1 1 /cv 1 1 2 1 1 /odd 1 1 2 
>2||/odd||l-^(^||/ev||l + 6||/odd||l) 

>~IIAvlli>~ll/lll 
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which shows that the lowest eigenvalue of A + P dd is no smaller than —1/18. 

□ 



According to Theorem [TJ then, if we choose potentials V Q dd and V ev which 
are equal to r~ 2 P odd , respectively r~ 2 P ev , for r > 1, and are equal to +1 for 
r < 1, then the first operator has a finite number of eigenvalues, while the 
second has an infinite number. On the other hand, it is clear that the region 
of phase space where a(L ev ) < —E is the disjoint union of two identical 
copies of the region where cx(L dd) < —E. So there are an infinite number 
of disjoint copies of the unit cube, under canonical transformations, in the 
second region if and only if there are in the first region. Moreover, this 
holds regardless of the precise conditions we place on these maps. Hence a 
Fefferman-Phong-type heuristic is not effective in estimating the number of 
eigenvalues of both these operators. 
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